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Abstract 

A general method for solving the so-called quantum inverse scattering prob- 
lem (namely the reconstruction of local quantum (field) operators in term of 
the quantum monodromy matrix satisfying a Yang-Baxter quadratic algebra 
governed by an i?-matrix) for a large class of lattice quantum integrable mod- 
els is given. The principal requirement being the initial condition (R(0) = P, 
the permutation operator) for the quantum i?-matrix solving the Yang-Baxter 
equation, it applies not only to most known integrable fundamental lattice 
models (such as Heisenberg spin chains) but also to lattice models with arbi- 
trary number of impurities and to the so-called fused lattice models (including 
integrable higher spin generalizations of Heisenberg chains). Our method is 
then applied to several important examples like the sl n XX Z model, the XYZ 
spin-Tj chain and also to the spin-s Heisenberg chains. 
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1 Introduction 



One of the main challenging problem in the theory of quantum integrable models, after di- 
agonalizing the corresponding Hamiltonians, is to find explicit and manageable expressions 
for their correlation functions. This is well known to be a very hard problem in general 
situations, even for quantum integrable models in two dimensions, and in fact besides the 
models which can be connected to free fermions or which are in the conformal field theory 
class, very few explicit solutions to this problem have been obtained so far. 

For example, in the case of the Heisenberg spin-^ chains Jjj], while the method for 
diagonalizing the Hamiltonian was designed in 1931 by H. Bethe || (for the isotropic XXX 
case), the first manageable expressions for its correlation functions (in the thermodynamic 
limit) have been given only in 1992 (in the massive regime) Q and then in 1996 (in the 
gapless regime) Q by Jimbo, Miwa and their collaborators, however based on some nice 
hypothesis, in particular, on the symmetry structure of the infinite chain. 

Motivated by this problem, it was shown for the first time in || that the explicit 
solution of the quantum inverse scattering problem (namely the reconstruction of any local 
spin operator at any site of the chain in terms of the elements of the quantum monodromy 
matrix satisfying a Yang-Baxter algebra, and containing creation/annihilation operators of 
Bethe eigenstates of the Hamiltonian) can be solved for this XXZ spin-^ model in a very 
simple (multiplicative) way. Moreover, using the algebraic Bethe ansatz framework [^,0], it 
was then used in J|] to derive explicit and simple expressions for the correlation functions of 
this model, even in the presence of a non-zero magnetic field. In the thermodynamic limit 
(the limit of infinite chain) it leads to multiple integral representations of these correlation 
functions. In the zero magnetic field limit, this approach gives then a complete proof of 
results and conjectures previously obtained in 

The elementary nature of the answer to the quantum inverse scattering problem, in 
this very representative example among the models solvable by means of the algebraic 
Bethe ansatz method, was quite unexpected. The algebraic Bethe ansatz method, also 
called quantum inverse scattering method, appeared 20 years ago as a quantum analogue 
of the classical inverse scattering problem approach to non-linear wave equations having 
soliton solutions, in order to solve quantum integrable models in two dimensions. The 
essential tools of this method are the quantum monodromy matrix satisfying quadratic 
commutation relations (Yang-Baxter algebra) which structures constants are given by an R- 
matrix solving the Yang-Baxter (cubic) equation |),1C]. It has been designed to diagonalize 
the corresponding Hamiltonians simultaneously with their associated commuting family of 
integrals of motion pLfTj. It is the quantum analogue of the direct part of the classical 
inverse scattering problem method in its Hamiltonian formulation [O], in which the Lax 
matrix is used to construct the monodromy matrix containing the action-angle variables 
which linearize the dynamics. However, the inverse scattering problem part of the classical 
theory, namely, the reconstruction of the local classical field variables contained in the 
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Lax matrix in terms of the elements of the monodromy matrix (and hence in terms of 
the action-angle variables) using the Gel'fand-Levitan-Marchenko equations [f[^-14|, being 



already a quite difficult problem to solve, it was not at all obvious to find a direct way to 
extend it to the quantum situation, although the motivations (form factors and correlation 



functions) were very clear from the very beginning [15-17]. 

The very remarkable feature of the solution to this problem given in |jj was not only 
that the quantum inverse scattering problem can be solved explicitely, but also that both 
its resulting expressions (reconstruction of the local spin operators at any site of the chain 
in terms of a simple product of the quantum monodromy matrix elements) and their proofs 
are very elementary. In turn, it essentially relies on the fact that the quantum i?-matrix 
i?(A, fi) solving the Yang-Baxter equation reduces to the permutation operator when the 
two spectral parameters A and \i are equal. 

This fact being almost a consequence of the Yang-Baxter equation itself, and hence 
satisfied for very generic cases, it immediatly suggested that the quantum inverse scattering 
problem can indeed be solved for almost all known quantum integrable (lattice) models^. 
This is to be compared to the tentatives of solving the corresponding problem directly for 
the continuum quantum integrable field theories such as the Sine-Gordon model, using a 
quantum version of the Gel'fand-Levitan-Marchenko equations, where it appeared to be 
an extremely difficult problem to handle, due in particular to the presence of boundstates 



lggOJ. 

The results we will present here show that for lattice models at least (whatever com- 
plicated is the corresponding spectrum) this problem admits a simple and purely algebraic 
solution. 

More precisely, let us consider a quantum lattice model associated to the finite dimen- 



sional representations of some quantum algebra A 22|, and let W r (X), where r + 1 is 
the dimension, such representations depending on a spectral parameter A (here we have in 
mind for example the finite dimensional evaluation representations of some quantum affine 
algebra). Let us suppose that the corresponding quantum i?-matrix R(\, fi) is invertible 
in any irreducible tensor product End(W rn (A) <8> W rn (fi)) and reduces to the permutation 
operator for equal spectral parameters. Then, let the quantum Lax operator L n {\) at site 
n defining this lattice model be obtained as some finite dimensional representations of the 
quantum i?-matrix in End(F(A) <g> W rn (£ n )) as, 

in (A - ^ + U ) = (A - U + U ) 

= R™ r »(\-£ n ) eEnd(V(\)®W rn (U)), 

where the auxiliary space V(A) (labelled by subscript 0) is for example a fondamental 
(evaluation) representation of A, W rn (£ n ) is the quantum local space at site n, and £ n is 

1 This idea was given by the authors in various conferences, including the"Seminaire de Mathematique 
superieure" held in Montreal in the summer 1999 |18|. 
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an inhomogeneity parameter attached to the site n, the space of states of the model being 
given by, 



H~W ri {Z 1 )®---®W rN (Zn)- 

Let L n (X) be choosen moreover such that we can associate to it an auxiliary quantum 
Lax operator C n {\) obtained by fusion 21, p3f-25] of L n {\) (meaning that its operator 
matrix elements are functions of the one's contained in L n {\)) and given by the quantum 
i?-matrix in End(W / rn (A) <S> W rn (<£ n )). Then, even in the general situation where the local 
quantum space of states W rn (£ n ) are different from one site to another, it is possible to 
give an explicit and elementary solution of the reconstruction of the local operators at any 
site n of the lattice in terms of elements of the quantum monodromy matrix associated to 
L n {\). In fact the resulting expressions for these local operators reconstructions and their 
proofs are very direct generalizations of the one's given for the Heisenberg XXZ case in 
which is quite generic in this respect. 

^From a more universal point of view, which will be useful when considering the exten- 
sion of our results to the field theory case for which infinite dimensional representations 
should be considered, the solution of the quantum inverse scattering problem relies in turn 
to the very simple fact that the associated universal 7£-matrix is the canonical (identity) 
element in the corresponding double Hopf algebra T>{A) |2q1 . 

The main purpose of this article is to give a simple and general solution to the quantum 
inverse scattering problem for any quantum integrable lattice model of the above type 
solvable by means of the algebraic Bethe ansatz method ||. This result applies to most 
known quantum lattice integrable models, even with arbitrary number of impurities at 
different lattice sites, since our method takes into account the possibility to have any 
chosen (finite dimensional) representations for the local space of states at each site which 
could be different for different sites. We will also give explicit formulas for several known 
interesting examples including not only the so-called fundamental lattice models (such as 
XYZ spin-^ or sl n XXZ models, for which the result is essentially a direct copy of the case 
described in ||), but also for the more general cases of fused lattice models (such as the 
higher spin generalisations of the Heisenberg chains). In that last interesting situation, we 
will also discuss the implications of these results for the classical inverse scattering problem. 

The application of these results to the actual computation of correlation functions of 
the above models will be considered in forthcoming publications. 

This article is organized as follows. In section 2, we first recall the results for the 
spin-^ XXZ Heisenberg spin chain as given in ||, which will be a guideline all along the 
article. Then the results are extended in a quite direct way to a large class of lattice 
quantum integrable models in section 3. In section 4, we apply these results to the so- 
called fundamental lattice models for which the quantum Lax operator is given essentially 
by the corresponding quantum i?-matrix, namely, for which in particular the auxiliary and 
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quantum spaces are isomorphic. In section 5, we consider the so-called fused lattice models, 
with the example of the XXX spin-s model. Some conclusions and comments are given in 
the last section. 



2 Solution of the quantum inverse problem for the Heisen- 
berg spin-1/2 chain 

In this section, we recall the solution of the quantum inverse problem which was given 
in H in the case of the XXZ spin-1/2 Heisenberg chain. 

The XXZ Heisenberg chain of length N corresponds to the following Hamiltonian, 



N 



rn+1 



(1) 



m=l 



where <r^, a = x,y,z, are the Pauli matrices acting in the two-dimensional space TL m at 
site m, A = cosh 77 is the anisotropy parameter, the particular case A = 1 corresponding 
to the XXX chain. We impose periodic boundary conditions. 

In the algebraic Bethe ansatz framework, the corresponding quantum i?-matrix is 



R(X,fx) 



where 



c(\,n) 



(I 0\ 

b(\,fi) c(A,/i) 

c(A,/i) b(\,fj,) 
\0 01/ 



(2) 



<p(v) 



with the function (p defined as 

99(A) = A in the XXX case, 

(f(\) = sinh(A) in the XXZ case. 



(3) 
(4) 



(5) 
(6) 



The i?-matrix is a linear operator in the tensor product of two two-dimensional linear spaces 
Vi (g) V2, where each Vi is isomorphic to C 2 , and depends generically on (the difference of) 
two spectral parameters Ai and A2 associated to these two vector spaces. It is denoted by 
^12 (Ai, A2) or by i?i2(Ai — A2). Such an i?-matrix satisfies the Yang-Baxter equation, 



i?i2(Ai,A 2 ) flis(Ai,A 3 ) i?2 3 (A 2 ,A 3 ) = i?23(A 2 ,A 3 ) i?i 3 (Ai,A 3 ) i?i 2 (Ai,A 2 ) 



(7) 
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Moreover, it is unitary, except for a finite number of values of the spectral parameter (such 
that 6(Ai,A 2 ) = ±c(Ai,A 2 )), 



i?i 2 (Ai-A 2 ) #2i(A 2 -Ai) = l, (8) 
and reduces to the permutation operator for a particular value of the spectral parameter, 

#12(0) =Pl2 

Identifying one of the two linear spaces in the .R-matrix with the quantum local space 
7i n at site n, one obtains the quantum L-operator of the model at site n as 

L n (\^ n ) = R 0n (\-( n ), (9) 

where £ n is an arbitrary (inhomogeneity) parameter dependent on the site n. The Hamilto- 
nian (|l|) corresponds to the homogeneous case where all the inhomogeneity parameters £ n 
are equals. The subscripts mean here that Rq u acts on the tensor product of the auxiliary 
space Vo — C 2 by the quantum space TC n at site n. The quantum monodromy matrix of 
the total chain is then defined as the ordered product of all L-operators along the chain, 

Ib(A) = r ,i..jv(A;£i,... ,£jv) = R 0N (\ - t; N ) . . . R 01 (\ - (10) 

and can be represented in the auxiliary space Vo as a 2 x 2 matrix, 

S3)- 

whose matrix elements A, B, C, D are linear operators on the quantum space of states of 

N 

the chain 7i = ® Ti n . Their commutation relations are given by the following relation on 
n=l 

C 2 ®C 2 : 

i?i 2 (A,M) Ti(A) T 2 (fi) = T 2 (n) Ti(A) i?i 2 (A,^), (12) 
with the usual tensor notations Ti(A) = T(A) <8> Id and T 2 (/i) = Id (g) T{fj). 

In the algebraic Bethe ansatz framework, the matrix elements of the monodromy ma- 
trix are used to construct the eigenstates of the Hamiltonian. In particular, the transfer 
matrices t(A), which are the trace -A(A) + -D(A) of the corresponding monodromy matrices 



(11), commute with each other for different values of the spectral parameter A, and the 
Hamiltonian (|l]) can be obtained from these transfer matrices (in the homogeneous case) 
by means of trace identities. In this context, B and C are used respectively as creation and 
annihilation operators of Bethe eigenstates. The common eigenstates of the Hamiltonian 
and the commuting transfer matrices are constructed by successive actions of operators 
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B(Xj) on a reference state | 0), which in this case is the ferromagnetic state with all spins 
up, for any set of n spectral parameters 1 ^ j ^ n} solution of Bethe equations. 

The quantum inverse problem consists in reconstructing the local operators, which here 
are the local spin operators = S x ±iS y , S z at a given site i of the chain, with S a = \o a , 
a = x, y, z, a a being the standard Pauli matrices, in terms of the matrix elements A, B, 
C, D of the monodromy matrix used to construct the eigenstates of the Hamiltonian. This 
was done in ||, where we obtained the following result : 

i— l n ( i n -l 



Si = \\\(A + D){u}B{i l )\i\{A + D)(i a )\ , (13) 

^ a=l ' a=l ' 

St = \ll^ + ma)}c(C l )\ll(A + D)^ a )\ \ (14) 

^ a=l ' a=l ' 

^={n^+^)(^)}^-^)(e*){n(^+^)(w} ( is ) 



In other terms, the elementary operators E n n ' n , e' n ,e n E {1,2}, acting on the local quan- 
tum space 7i n ~ C 2 at site n as the 2x2 matrix E\- ,£ = 5i £ '5j £ , 1 ^ i,j ^ 2 



arc 



reconstructed in the following way in terms of the matrix elements T £n ^ n of the quantum 



monodromy matrix (11): 



E ey n = I l[(A + D)(u}T £n>£ ^ k ) { f[(A + D)(^)\ \ (16) 

I. 01=1 ' 01=1 ' 

Note that the expressions (13) - (|l~6|), which are given here in the more general case 
with inhomogeneity parameters £j, are in particular valid in the homogeneous limit where 
all parameters £j are set equals to zero. 

These formulae were used in || to obtain an explicit determinant representation for 
the form factors of the finite chain, and in || to derive integral representations for general 
correlation functions of the chain in a magnetic field at the thermodynamic limit. 

The proof of (13)-(15) proposed in || was very simple. It relies essentially on the fact 
that these formulae can be easily proven, using that R(0) = P, for an operator acting on 
the first site on the chain (cf. Lemma 4.2 of ||), and that the problem can be reduced to 
this case by means of a propagator, which is a shift operator from site 1 to site i, and which 
is given as some product of i?-matrices or as a product of transfer matrices (cf. equations 
(4.3) and (4.4) of ||). In |j, the propagator was also expressed in terms of factorizing 
F-matrices, but it is clear that its properties rely actually only on its expression (equation 
(4.3) of H)) as a product of i?-matrices. Thus, the proof proposed in || is quite general, 
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and can be extended to a large class of lattice integrable models, as it will be shown in 
the next sections. In particular, its generalisation to fundamental models is quite obvious 
(see section 4 for some other interesting examples), but we will also show that one can 
obtain an explicit reconstruction for fused lattice models as well, such as for example the 
Heisenberg chains of higher spins (see section 5). 

3 General solution for lattice models 

In this section, we show that the resolution of the quantum inverse problem given in || in 
the case of the XXZ spin-1/2 Heisenberg chain is quite general and can be extended to a 
large class of lattice models. 

Let A be a Hopf algebra, for example a Yangian or a quantum affine algebra, which 
admits a family of finite dimensional evaluation representations W r (X), where r + 1 is the 
dimension, depending on a spectral parameter A. Suppose that, for any irreducible tensor 
product of two such representations W r (X) (8> W s (/jl), there exist a corresponding i?-matrix 
£>W r W s ^ _ ^ so i v i n g t ne Yang-Baxter equation, and which is unitary (except for a finite 
number of values of the spectral parameter). Suppose moreover that the matrix R WrWr (X), 
acting on the tensor product of any two spaces with same dimension r + 1, reduces to a 
permutation operator for a particular value of the spectral parameter : 

JjWrW r fQ\ _ pW r W, 

For notation conveniences, we will consider the additive spectral parameter case. However, 
it will be clear all along the results and proofs, that our method being purely algebraic, 
it extend to more general situations as well. Let us also suppose that there exists some 
fusion rules which allow the explicit construction, from the i2-matrices R (A) associated 
to fundamental representations V(X), of the i2-matrices associated to higher dimensional 
representations : 

R w r w 3{X) = T w T w s y [ R w^ _ (18) 

Then let us consider a quantum lattice model, solvable by means of algebraic Bethe 
ansatz, whose space of states is isomorphic to a tensor product of evaluation representa- 
tions: 

H~W ri (t 1 )®---®W rN (t n ). 

Here, W rn (£ n ) is the quantum local space at site n, and £ n is an inhomogeneity parameter 
attached to the site n. Let us define the quantum L-operator of the model at site n by 

L n (X -£ n + U ) = L^n (A — £ n + U ) 

= iCSA-£ n ) eEnd(V(X)®W rn (( n )), 
7 



(17) 



where vq € C, the auxiliary space V(A) (labelled by subscript 0) being a fondamental 
(evaluation) finite dimensional representation of A. The monodromy matrix, whose matrix 
elements are used to create eigenstates of the Hamiltonian in the framework of algebraic 
Bethe ansatz, is given by the ordered product of these L-operators along the chain: 

T 0) i..jv(A) = L 0N {\ -frf + vo)... L i(A - £1 + u ). (19) 

In this context, the quantum inverse problem consists in reconstructing local operators 
E%, 1 ^ a ^ (r n + l) 2 , which generate the local quantum space End(Ty rn ) at any given 
site n of the chain, in terms of the operator entries of the monodromy matrix (|lS|). 

The idea of this reconstruction, for an operator E° acting at a given site n of the chain, 
is to consider an auxiliary quantum C operator for which the dimension of the auxiliary 
space W is equal to the dimension of the local quantum space W Tn at site n, and such that 
this Lax operator at site n reduces to the permutation operator ( |17| ) in a particular value 
of the spectral parameter. Thus, we define the following auxiliary Lax operator for any site 
i of the chain, with an auxiliary matrix space W isomorphic to the space of states W rn , 
for a given n, 

4 n) (A-^ + ^o)=<^(A-^ 

acting in the tensor product W ® W n ~ W rn <8> W n , and such that it can be obtained 
by fusion from our starting quantum Lax operator L n . Here and in the following, the 
upperscript (r n ) means that the auxiliary space choosen to construct the corresponding 
object is isomorphic to W Tn . For this special choice of the auxiliary space, let us note 
T^ Vn \X) and t( r ")(A) the corresponding monodromy and transfer matrices: 

T^\X) = T$l N (\) = C$(\ -Zn + uq)... 4i n) (A - £i + uq), (20) 

t (r n ) {X) = t (r : ) N (A) = tr() (^"^(A)) . (21) 

We have the following proposition: 

Proposition 1. An operator E® G End(VF rn ) in a given site n of the chain can be ex- 
pressed in the following way in terms of the monodromy matrix T^?' N (£ n ) constructed 
with an auxiliary space W ~ W rn and of the transfer matrices constructed with the differ- 
ent auxiliary spaces isomorphic to the different quantum local space of states W Ti for any 
site index i, v- ri \^i), 1 ^ i ^ n: 

, n— 1 n (- n n — l 

K = \U tr (EgT^jtfa)) n^te) • ( 22 ) 

^ i=i ' < i=i ' 

where Eq denotes the operator E% acting in the auxiliary space W ~ W rn . 
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Proof. It comes directly from the two following lemmas. 

Lemma 1. The product of the above defined transfer matrices (which all commute among 
themselves due to the Yang-Baxter relation), 

n-l 

T^ n =n* (ri) te) 
i=i 

is an invertible shift operator from site 1 to site n, whose action on the monodromy matrix 
(for any choice of the auxiliary space labelled by 0' ) is given by 

,n...Nl...n-l 

(A). (23) 

Proof. The expression (p3|) can be proved by induction on n. The transfer matrix t^ n ^ N (^ n ) 
constructed with auxiliary space W isomorphic to W Tn can be expressed as a product of 
R-matrices: 

*i...iv(?n) = tr [R m N {U - £n) ■ ■ ■ R 0n+1 {U-U+i)P Q n 



tr ( 



X R nN &-W)-.-«„, n+ i (Cn-Cn+l)J 



w r „w riv , . TfryrnVYT n+1 , . 



where we used the cyclicity of the trace and the fact that, as W ~ Wr n j i?^^" (0) is the 
permutation P Qn rn of the auxiliary space W and the quantum local space at site n, W Tn . 
Using also the expression of the monodromy matrix, for an auxiliary space W of arbitrary 
dimension, as a product of i?-matrices, 



Try 



0',n..JVl...n- 



X .Kq,^ (A — §atJ . . . K , n (A - 



one can easily show, by using successively the Yang-Baxter equation 

RniHn — £i) Ro'i(^ — 6i) -Ro'n(A — £ n ) = -Ro'n(A — f n ) Ro'i{^ ~ £i) Rni(£,n ~ 
for i = n + 1, . . . , iV, 1, . . . , n — 1 that, 

4 r ™Ar(£«) To', n ...Afl...n-l(A) = Ib',n+l..JV'l...n(^) ^L™jv(^ n )' 

for any n, which achieves the proof. □ 
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Lemma 2 . If the auxiliary space W has same dimension as the quantum space W ri at site 
1, one has the following identity, for any operator E a £ End(W / ri ~ W): 

(n) (c ^\ _ ™ .(r-i) 



tr {E-T^ N (^))=E-t^ 1 ). 

Proof. Starting from the expression of as a product of i?-matrices, on uses the 

fact that i? 01 n (0) is the permutation operator Pq\, which acts in the following way on 

"01 -fy) — "01- □ 

It should be noted that we have considered here the general inhomogeneity case, but 
of course, our result applies and simplifies for the homogeneous case where all parameters 
£ n are equal. 

Proposition [l] gives thus the reconstruction of local operators at site n in terms of 
the matrix elements of the monodromy matrix the auxiliary space of which has the same 
dimension as the quantum space of states at site n. But the creation and annihilation 
operators used when solving the model by means of algebraic Bethe ansatz are elements 



of the monodromy matrix Ti...at(A) (19) the auxiliary space of which is a fundamental 
representation of the algebra A. One can use then the fusion rules ( |18| ) to reexpress the 
matrix elements of the monodromy matrix (20), which appear in the expression (^) of 
proposition |l| in terms of those of (19), solving the problem. We will see some example in 
the next sections. 



4 Fundamental lattice models 

In this section, we apply the general formula of section 3 in the special case of fundamental 
lattice models. Fundamental lattice models are characterized by the fact that the quantum 
space of states at any site n of the lattice is given by the fundamental representation of the 
corresponding quantum algebra and is isomorphic to the auxiliary matrix space. Hence, 
the Lax operator is just given by the quantum i?-matrix in the tensor product of the 
fundamental representation by itself and indeed reduces to the permutation operator in 
this tensor product when the spectral parameter is equal to zero, 

L n (\a n ) = LiY n ^-an) 

= RlY(\-£ n ) G End(V(X) ® V(£ n )), 
where V is a fundamental representation. The space of states is here 
H =i V(£i) (2) • • • (8) V(£ N ). 
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The quantum monodromy matrix is given by 

?b,i.„jv(A) = L 0n (X,( n ) . . . L i(A,£i), 
and the transfer matrix is defined as usual by its trace in auxiliary space Vo, 

*i..jv(A) = tr (T 0) i.„jv(A)) • 

As a consequence, the formula (22) of Proposition [l] applies directly to give the reconstruc- 
tion of any local operator E® acting at a given site n of the chain, 

K = {um}^o(EST 0tl ... N (U) {flm} > (24) 

where Eq denotes the operator E" acting in the auxiliary space. Note that the homoge- 
neous case is obtained as usual by putting all inhomogeneity parameters £ n to zero in the 
above formula. 



Let us stress here that it is a very direct copy of our original formula (16) given in || 
for the case of the XX Z spin-^ chain, which is completely generic in this respect. This 
formula is of course valid for both inhomogemeous and homogeneous models which are 
obtained simply in the limit where all inhomogeneity parameters £ n are set equal to zero. 

4.1 The sl n+1 XXZ model 

The quantum i?-matrix of this model, solving the Yang-Baxter equation and associated to 
a fundamental representation of U q {sl n+ i), is defined by the following (n + l) 2 x [n + l) 2 
matrix I27H291, 



R{n) = ejt e_H_ = 

e x q - e A q 1 

with, 

RW =Y^E aa ®E Pfi + qY^ E aa ® E aa + (q- q' 1 ) ^ E a/3 <g> E pa , 
R(") = ^ E aa ® Epf + q- 1 Y E aa ® E aa - (q - q~ l ) ^ E a/3 ® E^, 

where E a @ is the elementary (n + 1) x (n + l)-matrix with only one non-zero element, 
E"? = 5i >a 5j t /3. Note that, R(0) = P, the permutation operator. For homogeneous models, 
the quantum Lax operator being directly equal to this i?-matrix, the quantum monodromy 
matrix is given by, 

r ( l*r(A) = 45(A) . . . 4?(A)4?(A). (26) 
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The transfer matrices t^ n \X) = trolg" ^(A) commute among themselves for arbitrary 
values of the spectral parameters A, and one can define the local quantum Hamiltonian by 
means of trace identities, 

= (q-q- l )4rlnt^\X) ] -Niq + q- 1 ). (27) 
aX | a=o 

The simultaneous diagonalization of the transfer matrix and of the Hamiltonian is given by 



the so-called nested Bethe ansatz p0,31|. The solution of the quantum inverse scattering 



problem consists in reconstructing the local operators 

Ep p,Ep , l^e' p ,e p ^n + l, 

acting in the local quantum space at site p, TC P ~ C n , as the elementary matrices 

^ij ' £p = 3i,e' $j,e p i m terms of the elements of the quantum monodromy matrix N (X), 
whose elements Tij(A), 2 ^ j ^ n + 1 are used in particular to create Bethe eigenstates. 
We have directly from the above formula, 



E £ /' £p 



t^(0)Y 1 tr (E$' £ » 2*g >JV (0)) U n \0)\ \ (28) 



where the trace can be computed in terms of the matrix elements {T^ jv(0))e p ,ep °f the 
quantum monodromy matrix (26) to give finally, 



i#* = |tW(0)j P 1 (^?.jv(0)) |i (n) (0)) • (29) 

This expression (^) gives the reconstruction of the quantum local operators for the homo- 
geneous chain and is a direct generalization of (|l6|). It is clear from the general formulas 
of the preceding section that it is also possible to write a similar formula for the inhomo- 
geneous case as well. 

4.2 The XYZ spin-± model 

The XYZ model for spin-^ defined on a chain of lenght N is given by the following Hamil- 
tonian operator, 

N 

HxYZ = ^2 i J x a m a m+l + J y a m a m+1 + J z^m a m+l}- ( 30 ) 

771=1 
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(31) 



with 



b{u) 



sn(tt) 



(32) 



c{u) 



sn(u + rf) ' 
sn(fy) 



(33) 



d{u) 



sn(u + rf) ' 
A; sn(r/) sn(u) 



(34) 



where sn(ti) = sn(u; k) is the jacobi elliptic function of modulus k, ^ k ^ 1. The 
constants in the Hamiltonian are related to the parameters as, 



Then it is obvious that this i?-matrix reduces to the permutation operator P for u = 0. 
Moreover it is unitary up to a numerical factor that can be taken into account in all our 
computations. Defining the quantum Lax operator as being the i?-matrix itself, all the 
requirements for applying our above general formula are satisfied. Hence, we just obtain 
the reconstruction of any spin operator at site n to be given by the same formula as 
in the XXZ case where we have just to replace the operators A,B,C,D in ( |i~6|) by the 
corresponding one's from the XYZ model. 

5 Quantum inverse problem for higher spin XXX chains 

In this section, we consider the integrable generalization of the XXX model with local spin 
variables S" realizing a representation of s/2 of dimension / + 1 (spin 1/2). This model is 
given by a Lax operator Lq^(A), acting on the tensor product of a 2-dimensional auxiliary 
space Vq by the local quantum space Ji n ~ C /+1 at site n, 



J x = 1 + k sn 2 (rj) 
J y = 1 — k sn 2 (rf) 
J z = cn(r/) dn(r/). 



(35) 
(36) 
(37) 




(38) 
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As in the spin 1/2 case, the corresponding monodromy matrix, given by the ordered product 
of these L-operators with 2-dimensional auxiliary space along the chain, 



T ( ) 1 1 ) .. Ar (A) = L«(A)...4 1 1 ) (A), (39) 
can be written on the auxiliary space Vq as a 2 x 2 matrix, 
(i) (X) _(A(X) B(X)\ 

whose matrix elements are used to construct eigenstates of the Hamiltonian. The commu- 
tation relations ( |i~2|) of these matrix elements are given by the same i?-matrix (|2|) as in the 
spin 1/2 case. 

The Hamiltonian of this model is given by means of trace identities involving the 
transfer matrices i^(A) constructed from L-operators l£(A) with (Z + l)-dimensional 
auxiliary space. This Hamiltonian describes nearest-neighbour interactions, which follows 
from the fact that Lq2(0) reduces to the permutation Lbn- The L-operators Lq^(A) with 
higher dimensional (s + l)-auxiliary space can be constructed from (|38| ) by the fusion 



procedure [23|. More precisely, if P^„ M is the projector (symmetrizer) onto the (s + 1)- 
dimensional subspace Vr ai „ M \ of the tensor product V ai <g> ■ ■ ■ Cg V ag of s auxiliary spaces 
V ai ~ C 2 , one has 

^i.a s ),JA)=P+.. 03 4 1 1 ) ,n(A+^)... 

■ • • t%,n (A - 1±1 P 1 V) ■ ■ ■ 42n (A - S ~^V) PL.a s • (41) 

Remember that, in all this section, the dimension of the local quantum spaces is fixed to 
be I + 1, and that the upperscript (s) labels only the dimension of the auxiliary space (to 
be s + 1). 

The quantum inverse problem consists here in reconstructing the local spin operators 
in terms of the matrix elements A, B, C, D of the monodromy matrix (^) with 2- 
dimensional auxiliary space. From Proposition [l], we know how to reconstruct them in 
terms of the matrix elements of the monodromy matrix T^^(A) with auxiliary space of 
dimension Z + 1: 



s;; = {^(o)}"" 1 tro (ss4 l l. N (o)) {t w (o)}~ n . (42) 



The problem is thus to express the quantities 

AW(A)=tr (5 Q T ( i l 7V (A)) (43) 



14 



for any component a of the spin operator and for any integer s ^ I. To do this, one uses 

(s) 

the fact that the monodromy matrix Tq x n with (s + l)-auxiliary space can be obtained 

by fusion in terms of the monodromy matrix N with 2-dimensional auxiliary space, 
similarly as in (41): 



T (al...a s ),l...N^ ~ P ai...a a Ta 1 ]l...N{^+ 2 ^ 

■ • • r i'!i-Jv( A " £± ^") • ■ ■ (> - '-ity Pi...... («) 

This fusion procedure enables us to formulate the following lemma: 

Lemma 3. For any integer s ^ I, the quantity A« (A) defined by equation ( |43] ) has the 
following expression in terms of the elements (A) of the monodromy matrix ( f40|) and of 
the transfer matrices t^ s \\) corresponding to (s' + 1)- dimensional auxiliary space (s' sj: 



s - 1 



Ai s) («) = £ + A£> (u + ^=f^) t^ (u + 



fc=l (45) 

Proof. We will show the formula (^) by recursion on s, for s ^ I. First, note that it is 
satisfied for s = 1. Moreover, it can be easily shown also for s = 2. 

Suppose now that it is satisfied until a given s < I. When fusioning a spin s/2 with spin 
1/2 representation, one has the following block-triangular decomposition for the product 
of the corresponding monodromy matrices ]25l, [33]l : 



where the block structure is in accordance with the decomposition 



(46) 



v ai ® v a2 ~ v {aia2) e v {aia2) , y ai ~ c z , v a2 ~ 



2 \r r^t rps+l 

s+2 



"(aia 2 ) — ^ J "(aia 2 ) — ^ J 

and where the function fi(u) is the quantum determinant of the monodromy matrix (p0|). 
Taking the trace in (46), one obtains the well-known fusion equation for the transfer ma- 
trices: 

t {1) (u + i±%) &\u) = & +1 \u + |) + /,(„ + S -n) t(-D( u - |). (47) 
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In the same way, decomposing the component of the spin S% ia2 = + S" 2 into a x + 
^<oia2) on ^( a i a 2) ® ^{aia 2 )> one can take the trace of the quantity 

in (46), which leads to the following recursion relation for Aa\u): 

A« (« + ^17) * (s) (u) + (« + ^t/) A W («) 

= A(f +1 )(n + |) + /,(n + At,- 1 )^ - |). (48) 

One can then easily show, replacing Aa ^ (it) and A& ^ (u) by the corresponding expressions 
(|45| ) in (48), and using the recursion relation ( p7[ ) to reconstruct the transfer matrices which 
appear in each term of the sum, that one obtains also the expression ( (45|) for Aa +1 \u). □ 

Consequently, the solution of the quantum inverse problem for local spin operators is 
given by the following proposition: 

Proposition 2. The local spin operators S*, = S% ± iSn at site n have the following 
expressions in terms of the matrix elements of the monodromy matrix (]40| ) ; 

s- = (t^o)}" 1 

x {yj ( <")(|„) B (^I„) ( (-»(^I,)}{ t (') ( o,}-, (49) 

^ k = l ) 

s+ = {t««(o)} M 

x { £ f <'-"(|,) (^i,) } {<"»(0)}-' , (50) 



5f 



k=l 

i-1 



{(C)(0)} 



k fc=l (51) 

Note that the above formulae can be easily generalized to an inhomogeneous model. 

Remark 1. The transfer matrices t^ s '(u) can be expressed (using fl47|)) in terms of t^\u) 
and the quantum determinant (see (2^,^31 S3). 
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Remark 2. The proof given here being quite general and relying essentially on the recursion 
relation (46) with respect to the dimension of the auxiliary space, similar results can be 
obtained for XXZ and XYZ higher spin chains. 

It is quite interesting to examine the quasi-classical limit of the above formulas. This 
amounts to consider the limit in which the quantum parameter r\ goes to zero while the 
spin s of the quantum representation space goes to infinity, their product (the classical 
spin) being kept fixed. In the above formulas, the sum on spin and shifts becomes an 
integral in this limit. 

In this context, it is worth mentioning that while the solution to the quantum inverse 
scattering problem admits a very simple solution (essentially as a product of elements 
of the quantum monodromy matrix), the semi-classical limit of this expression (which 
is equivalent to the Gel'fand-Levitan-Marchenko equations) is more involved. This is in 
particular due to the fact that while every term in the multiplicative reconstruction formula 
for the quantum model becomes in this limit a function of the classical spin, and commutes 
with any other term in this expression, there is an "accumulation" of such terms which 
number becomes of the order s (i.e, of order - in the limit r\ — > 0). Hence, this gives non- 
trivial contribution (coming from non-commutative effects) in the classical limit , explaining 
in some sense why the classical case seems a priori much more involved than the quantum 
one. These effects results in effective "universal" Backlund transformations coming from 
the multiple dressing of the A, B, C, D operators in the quantum formulas by transfer 
matrices. It also shed some light on the reason why the direct quantization of the Gel'fand- 
Levitan-Marchenko equations is a priori a very difficult problem to deal with, since it would 
be very hard to guess the mutiplicative quantum formula from it. We expect do discuss 
this classical limit in a separate publication. 

6 Conclusion 

We have shown in this article that the quantum inverse scattering problem can be solved 
for almost all known quantum integrable lattice models with arbitrary finite dimensional 
local space of states (which could be different from one site to another). Our derivation 
relies on the unitaity of the quantum i?-matrix, its initial condition, R(0) = P, and on the 
existence of fusion rules relating quantum Lax operators associated to different auxiliary 
spaces. Hence, our results apply not only to the fundamental lattice models (such as the 
Heisenberg spin chains and their higher rank algebra generalisations) but also to lattice 
models with arbitrary number of impurities and to the so-called fused lattice models, like 
the higher spin integrable generalisations of Heisenberg chains. We would like to stress 
that besides some elementary technical details, the fundamental lattice model situation, 
both in the result and in its proof, is a very direct copy of our original derivation for the 
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XX Z spini case ||. For other more general models, fusion should be used, leading to nice 
dressed formulas such as (49)-(51). Another interesting case, that we have not considered 
here, concerns quantum integrable lattice field theory models, such as the Sine-Gordon 
model. In this context, we will have to deal with auxiliary quantum Lax operators having 
infinite dimensional representation auxiliary space, and also to define a trace operation 
for such representations. We expect the resulting reconstruction expressions for the local 
quantum fields to be given by formulas very similar to the infinite spin limit of the above 
result for spin-s Heisenberg chains. 

These results can now be applied to the computation of form factors and correlation 
functions of these models along the line we described in [|| . 

As a last remark, we note that after the completion of this manuscript, a paper, strongly 
inspired by our original work ||, and the seminar given by one of us Jl8|], appeared in the 



hep-th/preprint archive ( hep-th /9910253] ) , also dealing with the quantum inverse problem, 



however only in the particular case of fundamental graded lattice models. 
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